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ABSTRACT 

The innermost parts of dusty debris disks around main sequence stars are currently poorly known 
due to the high contrast and small angular separation with their parent stars. Using near-infrared 
interferometry, we aim to detect the signature of hot dust around the nearby A4V star Fomalhaut, 
which has already been suggested to harbor a warm dust population in addition to a cold dust ring 
located at about 140 AU. Archival data obtained with the VINCI instrument at the VLTI are used 
to study the fringe visibility of the Fomalhaut system at projected baseline lengths ranging from 4 m 
to 140 m in the K band. A significant visibility deficit is observed at short baselines with respect 
to the expected visibility of the sole stellar photosphere. This is interpreted as the signature of 
resolved circumstellar emission, producing a relative flux of 0.88% ± 0.12% with respect to the stellar 
photosphere. While our interferometric data cannot directly constrain the morphology of the excess 
emission source, complementary data from the literature allow us to discard an off-axis point-like 
object as the source of circumstellar emission. We argue that the thermal emission from hot dusty 
grains located within 6 AU from Fomalhaut is the most plausible explanation for the detected excess. 
Our study also provides a revised limb-darkened diameter for Fomalhaut (#ld = 2.223 ± 0.022 mas), 
taking into account the effect of the resolved circumstellar emission. 

Subject headings: Circumstellar matter — techniques: interferometric — stars: individual (Fomalhaut) 



1. INTRODUCTION 



The young (~200Myr, iDi Folco et ail 120041) and 
nearby (7.7 pc) A4 main sequence star Fomalhaut 
(a PsA, HD 216956) has been the focus of much at- 
tention during the last decade in the context of plane- 
tary system studies. Although the discovery of a cold 
debris disk around this bright star dates back to early 
observations wi th the InfraRed Astronomical Satellite 
(|Aumannl I1985D . the first resolved observations of its 
disk have been obt ained only in the lat e 90's in the sub- 
millimeter regime dHolland et alJll998ft. a nd a few years 
later at optical wavelengths (|Kalas et alj I2005T ). Fomal- 
haut has nowadays one of the best studied debris disk, 
with most of its (cold) dust arranged in a narrow ring 
at 140 AU from its host star. This ring shows interesting 
features, such as a sharp inner edge that has been in- 
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terpreted as the result of the gravitational influe nce of a 
mass ive planet located just inside the dust ring jQuillcn 
2006). The existence of this predicted planetary com- 
panion was recently confirmed b y direct observatio ns of 
Fomalhaut in the optical regime (jKalas et alJ l2008V fur- 
ther boosting the general interest of this system. The 
Fomalhaut debris disk has additionally been suggested to 
contain a population of warm dust thanks to partially re- 
solved Spi tzer/MIPS observation s and Spitzer/IRS spec- 
troscopy (|Stapelfeldt et al.ll200l . 

Although the Fomalhaut debris disk has been stud- 
ied in many details during the last few years, its inner 
dust content remains rather elusive, with only weak con- 
straints on a warm dust population within 20 AU pro- 
vided by Spitzer observations. The main challenges for 
characterizing this warm dust population are the small 
angular separation and the high contrast between the 
star and the inner disk. Infrared interferometry is an 
appropriate tool to tackle these challenges: by providing 
an angular resolution as good as a few milli-arcseconds 
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(mas), it can potentially resolve dust populations down 
to a fraction of an AU from Fomalhaut. This paper is the 
first of a series that aims at using infrared interferometry 
to characterize the warm dust content in the inner few 
AUs around Fomalhaut, and thereby provide a better 
view on the global architecture of its planetary system. 

The present paper focuses on the search for hot cir- 
cumstellar emission in the near-infrared regime with the 
VINCI instrument of the VLT Interferometer. The prin- 
ciple for warm dust detection with interferometry is 
based on the fact that the stellar photosphere and its 
surrounding dust disk have different spatial scales. For 
an A-type star at 7.7 pc, the angular diameter of the pho- 
tosphere is about 2 mas, while the circumstellar disk ex- 
tends beyond the sublimation radius of dust grains, typi- 
cally located around 0.15 AU (i.e., 20 mas) for black body 
grains sublimating at T su b — 1500 K. The circumstellar 
disk is thus fully resolved at baselines as short as 10 m in 
the near-infrared and contributes as an incoherent back- 
ground emission at longer baselines, while the stellar pho- 
tosphere is only resolved at about 150 m. The presence 
of the circumstellar disk then shows up as a decrease 
of visibility at baselines longer than about 10 m with re- 
spect to the expected visibility of the stellar photosphere. 
This detection method, which works best at short base- 
lines where the stellar photosphere is completely unre- 
solved, has already been successfully used to detect cir- 
cumstellar emissio n around the main sequence sta rs Vega 
(I Absil et alllgPOl . r Cet dDi Folco et all 12001^ Aql 
(jAbsil et all 120081 ) and (3 Leo (jAkeson et al.l l2009f i. In 
this paper, we revisit archival VLTI/VINCI observations 
of Fomalhaut and apply the same method to obtain first 
evidence for hot dust within 6AU from Fomalhaut. 

2. OBSERVATIONS AND DATA REDUCTION 

The bright star Fomalhaut has been observed on sev- 
eral occasions between 2001 and 2004 with VINCI, the 
VLT Interferometer Commissioning Instrument, which 
coherently combines the infrared light coming from two 
telescopes in th e infrared H and K b ands (for a detailed 
description, see lKervella et a l. 2003a). Some of these ob- 
servations have alrea dy be en reported in two paper s by 
iDi Folco et all (|2004[ ) and iLe Bouquin et alj (j2006f ). In 
this section, we describe these two data sets as well as 
an unpublished data set that we have extracted from the 
ESO archives. Most of the observations described here 
have been obtained with the 40-cm test siderostats of the 
VLTI. Their field-of-view (FOV) is limited by the use of 
single-mode fibers inside the VINCI instrument, and can 
be described by a 2D Gaussian function with a full width 
at half maximum (FWHM) of l'/6 in the K band under 
standard atmospheric conditions at Cerro Paranal. This 
translates into a linear FOV radius at half maximum of 
about 6 AU at the distance of Fomalhaut. 

2.1. Long-baseline data from \Di Folco et all 1200 A ) 

Observatio ns of Fomalhau t have been obtained in the 
K band by IDi Folco et alj (|2004l ) between November 
2001 and October 2002 with three baselines of vari- 
ous lengths and orientations: E0-G1 (66 m) and B3-M0 
(140m) with the 40cm test siderostats, and UT1-UT3 
(102m) with the 8m Unit Telescopes (see Fig. [T]). A 
few observations have also been obtained in the H band 
on the E0-G1 baseline, but because _ff-band data are 



not available on other baselines and therefore cannot be 
used to assess the presence of circumstellar emission in 
this particular waveband, these observations will not be 
further discussed here. 

In the present study, we will use the c a librat ed squared 
visibilities published bv IDi Folco et a l. (2004). Data re- 
duction was performed by using the wav elet analysis of 
the frin ge power spectrum described by iKervella et al.l 
(|2004aD and implemented in the VNDRS data reduc- 
tion software 2 . The resulting coherence factors (or raw 
squared visibilities) have been converted into calibrated 
visibilities by estimating the instrumental transfer func- 
tion with interleaved observations of calibrator stars with 
known diameters. The uncertainties on the resulting 
squared visibilities have been separated into statistical 
error (related to the dispersion of the raw visibilities) 
and systematic error (induced by the uncertainty on the 
calibrators' angular diamet ers). 

This data set has allowed IDi Folco et all (|2004D to de- 
rive an accurate estimation of the uniform-disk angular 
diameter (0 UD ) of Fomalhaut: 2.197±0.011±0.020mas. 
The first error bar refers to the statistical dispersion of 
the measurements whereas the second takes into account 
the systematic effects. This UD diameter can be con- 
verted in to an equivalent limb-darken ed diameter (#ld) 
following lHanburv Brown et al.l (|1974h : 

0ld(A) _ / l-ux/3 

0ud(A) V 1 - 7u x /15 ' [ ' 

with u\ the linear limb-darkening coefficient at a given 
wavelengt h A. U s ing th e tabulated limb-darkening coef- 
ficients of lClaretl (|2000f ). the resulting LD diameter reads 
2.228 ±0.011 ±0.020 mas. Due to the lack of short base- 
lines, IDi Folco et al.1 (|2004f ) could not thoroughly study 
the circumstellar environment of Fomalhaut, and were 
only able to derive an upper limit of about 3% on the flux 
ratio between the circumstellar and photospheric emis- 
sions. 

2.2. Medium-baseline da ta fro m \Le Bouauin et al\ 

Observations of Fomalhaut have been obtained with 
VINCI in the K band using two 40 cm siderostats on the 
64 m D0-H0 baseline (see Fig. [1}. The observations were 
performed on five consecutive nights from October 8th 
to October 12th, 2004. A large par t of this data set was 
reduced bv lLe Bouquin et"afl (|2006f ) using a custom data 
reduction procedure, and was subsequently used to val- 
idate the new integrated optics beam-combiner IONIC- 
2TK, which was installed in replacement of the origi- 
nal MONA fibered b eam-combiner in A u gust 2004. The 
study performed by ILe Bouquin et al.l (2006) shows a 
very nice consistency between the MONA and IONIC- 
2TK versions of the VINCI instrument, including very 
similar spectral transmissions. The study also demon- 
strates the validity of a Fourier-based estimator to reduce 
the IONIC-2TK data. 

These observations were not intended to study the cir- 
cumstellar environment of Fomalhaut, so that the ob- 
serving strategy was not optimal in that respect. In 

2 Publicly available on the ESO VLTI web s ite (see 
http : //www. eso . org/projects/vlti/instru/vinci/drs/ ). 
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Fig. 1. — Schematic view of the VLTI baselines used in this study. Data at long baselines (UT1—UT3 , B3 -M0, EO— Gl; in red) an d 
intermediate baselines (DO— HO, in blue) are taken from previous studies respectively by Di Folco et al. (2004) and Lc Bouquin ct al. (2006), 
while the data obtained at the shortest baselines (E0-G0, in green) are unpublished data from the ESO archives. 



particular, only one calibrator star (88 Aqr) was used 
throughout the five nights. This is expected to produce 
significant systematic errors on the calibrated visibilities, 
all the more that, with its estimated ang ular diameter of 
3.240 ± 0.057 mas (jRichichi et al.l l2009f I. this calibrator 
star is partly resolved on the D0-H0 baseline. However, 
because some of the observations have been obtained for 
projected baseline lengths as short as 20 m, this data set 
is still very useful to assess the presence of circumstellar 
emission. At such baselines, the stellar photosphere is 
almost unresolved (V 2 ~ 0.98), while the fully resolved 
circumstellar emission can induce a significant deficit of 
squared visibility. It must also be noted that, based on 
this large data set and ta king advantage of the hi gh in- 
ternal stability of VINCI, iLe Bouquin et alJ (|2006h have 
derived a UD angular diameter of 3.26 ± 0.06 mas for 
88 Aqr, which we will use in this study, as well as an 
independent estimation of the UD angular diameter of 
Fomalhaut (2 .23±0.07mas). in goo d agreement with the 
estimation of iDi Folco et all (|2004D . 

For the present study, we have retrieved the full origi- 
nal data set (5 nights) and reduced it with VNDRS 3.1, 
the latest version of the VINCI data reduction software. 
A total of 86 individual observations are available on 
Fomalhaut for the five nights. The output of the soft- 
ware consists in wavelet-based Fourier domain estima- 
tions of the squared coherence factors [x 2 for the target 
star (Fomalhaut) and its calibrator (88 Aqr). These "un- 
calibrated visibilities" are illustrated for a representative 
night in Fig.O Coherence factors must then be converted 



into squared visibilities (V 2 ) 
the relation: 

V 2 = 



for the scientific object by 
(2) 



T 2 



where T 2 is the interferometric transfer function (TF), 
i.e., the response of the system to a point source. The 
estimation of the TF is based on the interleaved obser- 
vations of the calibrator star, whose V 2 can be com- 
puted from the a priori knowledge of its angular diam- 
eter ($ud = 3.26 ± 0.06 mas). This V 2 must be com- 
puted for the wide spectral bandwidth on which the 
VINCI observations are performed, taking into account 
the actual spectrum of the star (88 Aqr, a K1III gi- 
ant) and the spectral transmission of the VINCI instru- 
ment. The spe c tral tra nsmission of VINCI is taken from 
iKervella et al.l (|2003bD . and we have used for 88 Aqr 
the tabulated i^-band spectrum for a Kl III star in the 
iPicklesI (|1998f ) stellar spectral flux library. The interfer- 
ometric TF is then estimated as T 2 = n 2 /V 2 , and is 
illustrated for one representative night in Fig. [2j The 
statistical and systematic error bars on this estimation 
are computed following Appendix C of IKervella et al.l 
(I2004al) . 

Once the interferometric TF has been estimated at the 
time of the calibrator observations, the value of the TF 
and its error bars must be estimated at the time of the 
actual scientific observations. We have used two differ- 
ent techniques to perform the interpolation of the TF 
between the calibrator data points: a linear interpola- 
tion between the two calibrator observations bracketing a 
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Fig. 2. — Squared coherence factors (/i 2 ) and interferometric 
transfer function estimations (T 2 ) for one representative observ- 
ing night on the D0-H0 baseline with VINCI-IONIC- 2TK. Filled 
symbols correspond to /x 2 data points for the scientific target (Fo- 
malhaut, red diamonds) and the calibrator target (88 Aqr, blue 
squares). The interferometric transfer function (TF) is represented 
by the empty blue squares, which are derived from the calibra- 
tor measurements taking into account its known diameter. The 
solid line is an interpolation of the TF, obtained by a third de- 
gree polynomial fit to its estimated values. The 1-cr statistical and 
systematic error bars on the estimation of the TF are respectively 
represented by dotted and dashed lines. 



given scientific observations (as in iKervella et ai1l2004af ). 
and a global polynomi al fit on the whole night (as in 
iLe Bouquin et al.l 12006). We have checked that the two 
methods give identical results within error bars, and de- 
cided to use the second one for the rest of the study. 
A third order polynomial was found appropriate to fol- 
low the (small) variations of the interferometric TF dur- 
ing the whole night, as the target and calibrator moved 
across the sky. With this method, the statistical error 
bar on the interferometric TF is computed locally at any 
given time by using the covariance matrix on the param- 
eters of the polynomial fit, which is available on output of 
the standard IDL routine poly_f it .pro. The systematic 
error bar, on the other hand, is computed globally with 
a weighted sum of the systematic error bars on all T 2 
estimations, taking into account the correlation between 
calibrators. In the present case, since only one calibra- 
tor was used throughout the night, the global systematic 
error is equal to the average of the systematic error bars 
on all T 2 estimations. 

2.3. Short-baseline data from the ESO archives 

Eve n though the proj e cted baseline lengths covered 
by the ILe Bouquin et aLl (|2006f ) data set are sufficiently 
short to study the circumstcllar environment of Fomal- 
haut, we have searched the ESO archives for calibrated 
VINCI observations on even shorter baselines. We have 
identified four interesting nights in 2003 (September 
27th, September 30th, November 2nd and November 
3rd), where observations of Fomalhaut interleaved with a 
calibrator star (88 Aqr) have been obtained on the 16 m 
E0-G0 baseline, using the 40 cm test siderostats and the 
MONA beam combiner. We have used VNDRS 3.1 to 
reduce this data set. One night (November 2nd) was 
discarded due to poor atmospheric conditions, and the 
remaining three nights produced 44 individual observa- 
tions of Fomalhaut of sufficient quality. 

Following the method described in Sect ion [2721 we have 



Fig. 3. — Same as Fig. |2]for one representative observing night 
on the E0-G0 baseline with VINCI-MONA. A second degree poly- 
nomial was used to fit the estimated values of the interferometric 
TF. 
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Fig. 4. — Sampling of the Fourier (u,v) plane for our complete 
data set, using the same color code as in Fig. [T] 



used the interleaved observations of the calibrator star 
88 Aqr to estimate the interferometric TF. In the present 
case, a second-degree polynomial was found convenient 
to model the TF, as illustrated for a representative night 
in Fig. [31 The scatter in the interferometric TF esti- 
mations in Fig. [3] may seem larger than in Fig. [2j but 
this is actually mostly due to the different scales for the 
two plots and to the smaller systematic error bars in the 
present data set. The reduced size of the systematic er- 
ror bars is related to the fact that the calibrator star is 
mostly unresolved on the very short E0-G0 baseline, so 
that the uncertainty on its angular diameter translates 
into a much smaller systematic error on the TF estima- 
tion. 

Our final VINCI data set on Fomalhaut, spanning a 
large range of baseline lengths and azimuths, comprises 
a total of 171 calibrated squared visibilities. The spatial 
frequencies sampled by our observations are well spread 
in the Fourier (u, v) plane, as illustrated in Fig. 2J 

3. EVALUATING THE AMOUNT OF CIRCUMSTELLAR 
EMISSION 
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Fig. 5. — Result of the fit of an oblate limb-darkened stellar photosphere to our full data set. For the sake of clarity, the solid line 
corresponds to a circular photosphere with a diameter equal to the geometric mean of the diameter of the best-fit oblate stellar model, while 
the residuals are computed with the full 2D photospheric model. The dashed and dotted lines in the bottom panel represent respectively 
a la and a 3<j deviation with respect to the best-fit model. 



Let us first assess whether a realistic stellar photo- 
spheric model for Fomalhaut could reproduce the whole 
data set. The description of our oblate limb-darkened 
photospheric model is given in Appendix^ showing that 
only one parameter needs to be fitted (#ld, the geomet- 
ric mean of the limb-darkened diameter). Fitting our 
model to the whole data set gives 6*ld = 2.292 ± 0.013 ± 
0.050 mas (Fig. [5]), where the first error bar accounts for 
the statistical dispersion of the data while the second 
error bar is related to the systematic error in the eval- 
uation of the interferometric TF . This result is withi n 
2a of the previous estimation by iDi Folco et al.l {2004), 
but the residuals of the fit in Fig. [5] and the large re- 
duced chi square (\r = 2.43) clearly demonstrate that 
an oblate limb-darkened photosphere cannot reproduce 
the data satisfactorily. It must be noted that our data 
sample various baseline azimuths, which makes it diffi- 
cult to represent the model together with the data on a 
single plot. The solid line of Fig. [5] is actually computed 
for a perfectly circular photosphere of diameter #ld and 
is only an approximation of the actual model used to 
compute the fit residuals, which take into account the 
various baseline azimuths and the stellar oblateness. 



Accepting the photospheric posit ion angle recently 
measured bv lLe Bouauin etlrtl (|2009h . possible improve- 
ments of the fit could be obtained by changing either 
the limb-darkening coefficient (uk) or the photospheric 
oblateness (p) . This is investigated in Table [TJ which 
gives the goodness-of-fit for various values of these two 
parameters. The fit is slightly better for extreme limb- 
darkening coefficients (uk ^1), and for a slightly more 
oblate photosphere than expected from our model. But 
the differences in \r are so small that we can safely state 
that our interferometric data do not constrain the shape 
of the photosphere in a significant way. We will therefore 
use the most plausible values (oblateness p = 1.021 and 
position angle PA = 156°, see Appendix |A|) in the rest 
of this study. 

3.1. Fitting a star-disk model to the complete data set 

The failure of a realistic stellar photospheric model to 
reproduce the interferometric data urges us to increase 
the complexity of our model. The apparent decrement 
of visibility observed at short baselines and the slope in 
the residuals of the fit suggest that another source of 
emission, fully resolved by the interferometer, must be 



6 



Absil et al. 



TABLE 1 

GOODNESS-OF-FIT (Xr) F0R THE PHOTOSPHERIC MODEL AS A 
FUNCTION OF OBLATENESS AND LIMB-DARKENING 





u K = 0.0 


u K = 0.19 a 


«k = 0.5 


u K = 1.0 


p = 1.00 


2.46 


2.46 


2.45 


2.43 


p = 1.021 b 


2.44 


2.43 


2.42 


2.41 


p = 1.05 


2.42 


2.42 


2.41 


2.40 


p = 1.10 


2.46 


2.46 


2.45 


2.44 



a Expect ed value of t he limb-darkening coefficient % in the 
K band l lClaretll2000h . 

b Expected value of the photospheric oblateness p based on 
our oblate photosphere model (see Appendix [A} . 



present within the interferometric field-of-view. As a 
first step, we assume that this resolved emission is as- 
sociated with the circumstellar debris disk, and to repre- 
sent its contribution we use a simple model of a diffuse 
source uniformly distr ibuted acro s s the whole field-of- 
view. As discussed by Absil et al.l ([2006), such a model 
is a good approximation provided that the circumstellar 
disk is fully resolved at all baselines. It must be noted 
that in the present case, due to the very short baselines, 
this condition might not be completely fulfilled. The uni- 
form emission model should nevertheless provide a good 
estimation of the flux ratio between the integrated cir- 
cumstellar emission and the stellar photosphere. More 
realistic circumstellar emission models will be discussed 
in Section [5l 

In Fig. [SI we have fitted the whole VINCI data set with 
a limb-darkened oblate stellar photosphere surrounded 
by a uniform circumstellar emission. The quality of the 
fit is very satisfactory, with \r — 0-95 and no obvi- 
ous trend in the residuals, which are now nicely spread 
around 0. The best-fit mean limb-darkened diameter 
is f?LD = 2.218 ± 0.009 ± 0.042 mas, while the best-fit 
flux ratio between the circumstellar disk and the star is 
0.88±0.05±0.11%. Assuming that our model is represen- 
tative of the actual brightness distribution, circumstellar 
emission is thus detected at a 7a level. Another repre- 
sentation of our detection is given in Fig. [7] and discussed 
in Appendix [B] 

Systematic errors dominate the noise budget in our 
detection. This is not surprising as our short- and 
intermediate-baseline data are calibrated with a single 
calibrator star (88 Aqr). The uncertainty on this cal- 
ibrator's diameter is the main contributor to the sys- 
tematic error budget. Using a few different calibrators 
would have been preferable to reduce systematic errors, 
but the archival data were not originally meant to be 
used for this particular purpose. We are nonetheless con- 
fident in the robustness of our result for three main rea- 
sons. First, 88 Aqr is a K giant with a slow rotational 
velocity ( 3 . 6 km s ~ 1 . iHekker k. Melende d [2007T) that has 
already been used as an infrared photometric standard 
(jBouchet et al.lll991l ) and as an interferometric calibra- 
tor by various authors without any hint of unexpectedly 
high or low visibilities. This includes a recent study 
dedica ted to interferometric calibrators bv lRichichi et al.l 
(2009). where a UD diameter of 3.240 ± 0.057mas is 
derived. Second, we have independently estimated its 
angular diameter usin g surface-brightness relationships 
(jKervella et al.ll2004bh . giving a UD diameter of 3.31 ± 
0.08 mas in good agreement with the value used in this 



study (3.26 ± 0.06 mas). Finally, we note that an im- 
perfect model for 88 Aqr cannot be at the origin of the 
observed visibility deficit for the following reasons. Ar- 
tificially decreasing its angular diameter would increase 
the calibrated visibilities of Fomalhaut. This would how- 
ever not reconcile the Fomalhaut calibrated data set with 
a simple photospheric model, as it would have a much 
larger influence on the data collected at mid-range base- 
lines (D0-H0) than on the short baseline data (E0-G0). 
A change in the limb-darkening coefficient would not help 
either. The presence of a circumstellar environment (disk 
or companion) around the calibrator could only reduce 
the measured calibrator visibilities with respect to their 
expected values, which would result in an increase of the 
calibrated visibilities for the scientific target and thereby 
reduce the measured disk/star flux ratio with respect to 
its actual value. Therefore, the probability that our de- 
tection is based on a systematic effect related to the cal- 
ibrator is considered to be extremely low. 

Our revised estimation of the stellar limb-darkened di- 
ameter of Fomalhaut (0 LD = 2.218 ± 0.009 ± 0.042 mas) 
is within the e rror b ar of the previous estimation by 
iDi Folco et alJ (|2004l ). as it is mainly based on long- 
baseline data which are hardly affected by the presence of 
the circumstellar emission. The large systematic error is 
due to the fact that most data (at short- and mid-baseline 
lengths) have been calibrated with the same reference 
star. A better accuracy on the stellar diameter can actu- 
ally be obtained by fixing the disk/star contrast at 0.88% 
and by fitting only the long-baseline data, which have 
been obtained with five different calibrators. The final 
result is then 6»ld = 2.223±0.006±0.021mas (\l = 0.31), 
w hich is in very good agreement with the estimation 
of IDi Folco et all |2004). The actual angular diameter 
ranges between 2.246 ± 0.022 mas and 2.200 ± 0.022 mas 
(apparent major and minor axes), taking into account 
the apparent oblateness of 1.021. 

3.2. Fitting a binary star model to the complete data set 

Besides a circumstellar disk, another potential source 
of visibility deficit at short baselines would be a faint 
point-like object within the interferometric FOV around 
the target star. To reproduce the observed visibility 
drop, the off-axis object should have a flux ratio of about 
0.88% with respect to Fomalhaut as seen through the 
interferometer in the K band, which is equivalent to a 
magnitude K = 6.1 for a close companion. Because of 
the Gaussian shape of the off-axis coupling efficiency into 
the single-mode fibers, the companion must be brighter 
than K = 6.1 to reproduce the VINCI near- infrared ex- 
cess if its angular separation is a significant fraction of (or 
larger than) the radius at half maximum (0'/8 ) of the in- 
strum ental FOV. Following the discussion in Absil et al.l 
( 2006), we estimate that the presence of such a bright 
background object within the instrumental FOV is very 
unlikely (probability ~ 10~ 6 ). Therefore we only con- 
sider the case of a bound low-mass companion. 

To fit our interferometric data set, we construct a 
model of a binary star, assuming that the companion 
orbits within the plane defined by the outer dust ring 
imaged at visible and sub-millimetric wavelengths. We 
therefore set the inclination and position angle of the 
projected orbit to 65?9 and 156°, respectively. In prac- 
tice, coplanarity with the outer dust ring is probably not 
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Fig. 6. — Result of the fit of a star-disk model to our full data set. The solid line represents the best fit star-disk model, while the dotted 
line represents the best-fit result with a single star for comparison, using the same conventions as in Fig. [5] 



mandatory to ensure its stability, but the purpose here 
is only to check whether binary star models could actu- 
ally fit the data rather than to explore all the possible 
orbital solutions. For the sake of simplicity, we will fur- 
ther assume that the orbit is circular, so that we are left 
with three parameters: the semi-major axis, the orbital 
phase at a given time to and the binary flux ratio (or 
contrast). Using a large number of potential values for 
these three parameters, we compute the position of the 
companion at our observing dates and the associated vis- 
ibility of the binary system. We deduce the chi square 
between the observed visibilities and the computed ones, 
and for each couple semi-major axis / contrast, we search 
for the minimum \ 2 as a function of the orbital phase at 
to- Thereby, we produce a \ 2 map, which is represented 
in Fig. [5] as a function of semi-major axis and binary 
contrast. 

The x 2 map shows that a whole range of semi-major 
axes and contrasts fit the interferometric data in a sat- 
isfactory way. The effect of the Gaussian beam pro- 
file is evident in this figure, as the binary contrast re- 
quired to fit the data increases with the binary sepa- 
ration. We have restricted the plot range to contrasts 
smaller than 5%, as brighter companions would most 
probably have been detected by simple near-infrared 



spectro-photometric measurements. Because the instru- 
mental transmission becomes extremely small at linear 
radii larger than 30 AU (i.e., angular separations larger 
than 4"), binaries with larger semi-major axes cannot 
reproduce our interferometric data set. Conversely, for 
any semi-major axis smaller than about 25 AU, one can 
find suitable combinations of orbital parameters and flux 
ratio to fit the data satisfactorily (reduced x 2 ~ !)• 
For the closest companions, located within the radius at 
half maximum of the FOV transmission (ft' 8, i.e., about 
6AU), the best-fit binary flux ratio is about 0.88% as 
expected. 

In conclusion, even though our interferometric mea- 
surements span a large range of time, baseline lengths 
and azimuths, they are not sufficient to discriminate be- 
tween a circumstellar disk and a point-like companion as 
the source of visibility deficit. This is mostly due to the 
(lack of) observing strategy, which was not optimized for 
such a goal. 

4. NATURE OF THE NEAR-INFRARED EXCESS SOURCE 

The previous section has shown that our interferomet- 
ric data set is not sufficient to discriminate between a 
point-like object and an extended source as the origin of 
the near-infrared excess emission. Therefore, we need to 
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Fig. 7. — Result of the fit of a star-disk model to our full data set, showing the data for two particular baselines as a function of the hour 
angle of the observation (left: the small E0— GO baseline, right: the mid-range DO— HO baseline). The blue curve represents the visibility 
of the star alone for our best-fit star-disk model, and its thickness corresponds to the 3a uncertainty related to the final error bar on its 
diameter. The black solid line represents the best-fit star-disk model, with its 3<r confidence interval shown as dotted lines. 
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Fig. 8. — Reduced chi square map for the fit of a binary star model, plotted as a function of the orbital semi-major axis of the companion 
and of the binary flux ratio. The orbital dynamics are taken into account in this model, assuming a circular orbit, and the initial orbital 
phase is optimized to minimize the Xr- We have limited the scale range to Xr — 4 for the sake of clarity. Contours are shown for Xr 
between 1 and 3.5, with a step of 0.5. 

use complementary data to further constrain the nature 
of the near- infrared emission around Fomalhaut. In this 
section, we focus on the possible presence of a low-mass 
companion, as well as on potential non-thermal sources 
of extended emission. The thermal emission from a cir- 
cumstellar dust disk will be briefly discussed in Sections [5] 
and[6l 



4.1. Further constraints on a low-mass companion 

Since a low-mass companion cannot be ruled out by 
the interferometric data, we have searched the literature 



for possible additional constraints on its existence. Three 
main type of observations can be used in this context: ra- 
dial velocities, astrometric measurements, and high dy- 
namic range single-pupil images. Before discussing the 
existing data in the literature, it is useful to convert 
the companion flux ratio i nto an estimated ma ss, using 
the evolutionary models of iBaraffe et alJ (|1998f ) for low- 
mass stars. Assumin g an age of 200 Myr for Fomalhaut 
(|Di Folco et al.ll2004j) . and taking into account its dis- 
tance of 7.7 pc, a if -band magnitude of 6.1 for a close-in 
low-mass companion translates into a mass of 0.35 M@. 
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Wider companions would have larger masses to account 
for the reduced off-axis transmission of VLTI/ VINCI, 
e.g., 0.48 Mq at fX'8 (i.e., 6 AU) where the interferometric 
transmission has dropped by a factor 2. This typically 
corresponds to spectral types in the range M1V-M3V. 

Radial velocity monitori ng of Fomalhaut was r ecently 
performed with HARPS bv lLagrange et ail (|2009f ). in the 
context of a radial velocity survey of early-type main 
sequence stars. Their observations of Fomalhaut span 
a 325 day period. Although Fomalhaut was found to 
be variable in radial velocity, no companion was de- 
tected with the following 68.2% probability upper lim- 
its on companion masses as a function of orbital period: 
3.6Mj up at 3 days, 3.2Afj up at 10 days, and 10Mj up at 
100 days, the latter corresponding to a semi-major axis 
of 0.5 AU. Since such planetary companions would pro- 
duce flux ratios largely below 0.88% in the K band, this 
basically rules out any companion closer than 0.5 AU as 
the source of the detected near-infrared excess. 

An anomalous proper acceleration of 6.6mas/yr 2 has 
been obse rved in HiPPARCOS astrometry towards Fo- 
malhaut (jChiang et al.l I2009T ). This acceleration is 
marginally significant (~ 2a) , the Fomalhaut astrometry 
being otherwise stable. Such quasi-steady acceleration 
could be caused by a companion whose orbital period is 
longer than ~ 3 yr (i.e., semi-major axis a > 2.6AU).The 
required mass of this companion to explain the observed 
acceleration would be about 9 Mj up at 2.6 AU, and would 
increase to 50M Jup at 7 AU, or 300M Jup (i.e., 0.29 M ) 
at 15 AU. The latter would produce a if-band flux ratio 
of about 0.6% with respect to Fomalhaut. Companions of 
larger masses at the same orbital distances would be in- 
compatible with the HiPPARCOS data, while companions 
with orbital periods smaller than 3 yr should have masses 
smaller than 9 Mj up to be consistent with HiPPARCOS ob- 
servations. Because a companion with a 0.6% flux ratio 
at an orbital semi-major axis of 15 AU is mostly incon- 
sistent with the VINCI data (xl ~ L8 ' see Fi S- EJ) ? these 
astrometric measurements rule out the low-mass binary 
scenario as an explanation for the VINCI near-infrared 
excess up to orbital semi-major axes of about 15 AU. At 
larger distances, the mass of a bound companion com- 
patible with the HiPPARCOS acceleration would become 
marginally consistent with the VINCI near infrared ex- 
cess. 

High contrast imaging with large diffraction-limited 
telescopes can be used put complementary constraints 
on low-mass companions at larger angular distances. In 
particular, the HST/AC S observations of Fomalhaut by 
IKalas et al.l (|2005l l2008h show no low -mass companions 
at the sensitivity limit of the instrument, except for the 
planetary-mass companion at 119 AU (well o utside the 
VINCI field). According to IKalas et all (12001) . the ACS 
instrument can detect magnitude differences up to 11 at 
0'.'9 (size of the Lyot occulting spot) for the F606W filter. 
Assuming a similar dynamic range in the F814W band, 
this gives a lower limit / = 12.0 for the absolute mag- 
nitude of the companion in the / band, which roughly 
corresponds t o a 0.08 object at 200 Myr and 7.7 pc 
following the Baraffe et al. (1998) models. This object 
would be much too faint to explain the X-band excess 
derived from VINCI observations, so that the effect of 
a companion located at an angular distance larger than 
0'.'9 can be ruled out. For a low eccentricity orbit seen 



face-on, this corresponds to an orbital semi-major axis of 
7 AU. However, assuming the orbital plane of the poten- 
tial companion to be coplanar with the circumstellar disk 
(65?9 inclination) and taking into account the projection 
effect, the planet could be on a circular orbit with a semi- 
major axis as large as 14.3 AU while still being hidden 
under the ACS coronagraphic mask for a 1.73 yr time 
span, which co rresponds to the interval between the two 
ACS images of IKalas et all (f2005L [2008h . The HST/ACS 
images can thus rule out the binary scenario as the expla- 
nation of the VINCI near-infrared excess for semi-major 
axes of 14.3 AU and above. 

More recently, M band observations with the AO- 
assisted Clio imager on the 6.5-m MMT have further 
constrained the range of allowed com panion masses in 
the 5 -40 AU region around Fomalhaut (|Kenworthv et all 
2009). In particular, objects with masses larger than 
13Mj up can be ruled out at distances between 8 and 
40 AU. Taking into account all these constraints, the 
presence of a bound low-mass companion around Foma- 
lhaut, whatever its angular separation, cannot be at the 
origin of the near-infrared excess detected with VINCI. 
And since the probability to find a K ~ 6 background 
object in the very small VINCI field is negligible, we are 
left with an extended source as the only plausible origin 
for the detected near-infrared excess. 

4.2. Possible sources of extended emission 

Among the possible sources of extended emission in the 
close vicinity of Fomalhaut (mostly within 0'.'8), hot dust 
grains producing thermal emission in the near-infrared is 
the most straightforward scenario, since Fomalhaut is al- 
ready known to be surrounded by large quantities of cold 
dust as well as by a wa rm dust population within 20 AU 
()Stapelfeldt et alj |2004). However, the thermal emission 
from hot dust grains is not the only possible source of 
extended near-infrared emission around Fomalhaut. In 
particular, non-thermal emission could be produced by 
free-free radiation from an ionized wind or a gaseous disk. 

The possible presence of intense stellar winds around 
A-type mai n sequence stars ha s already been thoroughly 
discussed in lAbsil et alj (|2008l ). showing that such stars 
have generally weak winds that are not expected to pro- 
duce significant near-infrared excesses. The a bsence of 
detec table X-ray emission towards Fomalhaut (|Schmittl 
|1997|) further suggests that a strong stellar wind cannot 
be at the origin of the observed if-band excess. Gaseous 
disks, on the other hand, are generally occurring in classi- 
cal Be or Ae stars, a category of hot rapidly rotating stars 
close to the main sequence that exhibit prominent hydro- 
gen emission lines. Although Fomalhaut does not rotate 
as fast as classical Ae/Be stars, and has a spectral type 
later than most emission-line stars, the possible presence 
of a gaseous disk cannot be ruled out a priori. Since 
infrared excesses in such stars are generally correlated 
to the equivalent width of t heir Ha emission line (see 
e.g. Ivan Kerkwiik et al.|[T995L and references therein), we 
have searched the literature for possible evidence of Ha 
emission around Fomalhaut. The Ha photospheric ab- 
sor ption line profil e of F omalhau t was inspect e d in d etail 
bv [Gardiner et alj (|1999fl and ISmallev et all (|2002l ). In 
particular, the observed Ha line profile was fitted with a 
model Balmer profile to deri ve the effective stella r tem- 
perature, which is found bv ISmallev et alj (|2002l ) to be 
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in agreement with the fundamental effective tempera- 
ture derived by other methods. No peculiar feature is 
reported in the Ha absorption line, including no obvious 
emission line inside it. Furthermore, the photometric in- 
de x of the Ha line was moni t ored b etween 1968 and 1974 
bv lDachs fc Schmidt-Kalerl (|1975j ). showing no variabil- 
ity and no evidence of Ha emission for Fomalhaut in 19 
observations. 

Finally, let us note that Fomalhaut is located well 
inside the Local Interstellar Bubble (LIB), so that 
the heating of interstellar material proposed to explain 
the presence of infrared excesses around sh ell stars or 
A Bootis stars located abou t 100 pc away ([Abll200l 
iMartmez-Galarza et al.ll2009f ) is unlikely to be at the ori- 
gin of the observed near-infrared excess. Consequently, 
we consider the presence of hot dust in the inner plan- 
etary system of Fomalhaut as the only plausible expla- 
nation to the observed near-infrared excess, and further 
discuss this scenario in the next sections. 

5. CONSTRAINING THE INNER DISK MORPHOLOGY 

We have previously derived in Section 13.11 the inte- 
grated flux ratio between the circumstellar disk and the 
stellar photosphere using a very simple brightness distri- 
bution for the disk: a uniform emission across the instru- 
mental FOV. Accepting that an inner dust disk is indeed 
at the origin of the detected near-infrared excess, we dis- 
cuss here the limitations to this preliminary analysis and 
investigate whether the morphology of the disk could be 
constrained based solely on our interferometric data set. 

As a first step, we have tried to use three very differ- 
ent morphologies for the circumstellar disk to reproduce 
the data set: a uniform emission extending across the 
whole field-of-view as in Section 13.11 a narrow ring of 
dust 3 located at the sublimation radius (r su b — 0.1 AU 
for a sublimation tem perature T sl ,h = 170 K), and the 
zodiacal disk model of iKelsall et alj (|1998f ). which is im- 
plemented in the Zodipic package 4 . Fig. [9] illustrates 
the influence of the disk morphology on the squared vis- 
ibility at short baselines. The small angular extent of 
the ring structure leads to significant oscillation in the 
squared visibility, whose frequency depends on the az- 
imuth of the baseline since the disk is inclined. Smaller 
oscillations are visible at a similar frequency for the zo- 
diacal disk model, because the inner rim dominates the 
-fT-band emission of the disk. Fitting these three models 
to our complete data set leads to similar results in terms 
of best-fit flux ratio and reduced x 2 ( see Tabled), except 
for the ring structure which gives a larger best-fit con- 
trast of ~ 1.3% with a slightly increased x 2 - The larger 
best-fit contrast is mainly due to the fact that the ring 
is not fully resolved at the shortest baselines (< 10m), 
so that a larger disk emission is required to produce the 
same squared visibility deficit at those baselines. Note 
however that the best-fit contrast is only about 3tr above 
our original estimation, which indicates that, even with 
this extreme morphology, the final disk/star flux ratio 
does not heavily depend on the distribution of the cir- 

3 The ring width is forced to be much smaller than the interfer- 
ometer's angular resolution so that it does not impact the visibility 
estimations. 

4 Zodipic is an IDL program developed by M. Kuchner et al. for 
synthesizing images of exozodiacal clouds. It can be downloaded 
at http : //asd . gsf c . nasa . gov/Mar c . Kuchner /home . html 
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Fig. 9. — Comparison of our data set with various disk models at 
short baselines. The data have been separated into three ranges of 
baseline azimuths to allow the comparison, since the model visibil- 
ities depend on both the baseline length and azimuth. In our mod- 
els, the narrow ring is located at the sublimation radius of black 
body dust grains, assuming a subl imation temperatur e of 1700 K, 
while the zodiacal disk follows the Kclsall et al. (1998) model. 



TABLE 2 

Best-fit flux ratio and goodness of fit for three 
different disc models 





Uniform emission 


Zodiacal disc 


Narrow ring a 


Flux ratio 


0.88% 


0.85% 


1.26% 


*5 


0.95 


0.93 


1.24 



a The ring is located at the dust sublimation radius (~ 0.1 AU). 



cumstellar emission. 

As a second step, we have generated images of a geo- 
metrically and optically thin debris disk toy model based 
on pure black body assumption (grain temperature pro- 
portional to r~ a5 , with T = 1700K at 0.1 AU), using a 
range of values for the two most important geometric pa- 
rameters: the inner disk radius r- m and the exponent a of 
the power-law describing the density decrease as a func- 
tion of distance (n(r) oc r~ a ). The synthetic images only 
include thermal emission, which is expected to be largely 
dominant for the dust temperatures explored here. For 
each couple of parameters (r; n , a), we compute the vis- 
ibility of the star-disk system at the relevant observing 
dates and adjust the disk/star flux ratio in order to min- 
imize the x 2 distance between the modeled and observed 
visibilities. By exploring a range of possible values for 
the two parameters, we produce a x 2 map, which is dis- 
played in Fig. [10] together with the associated best-fit 
flux ratios. In these toy models, we have assumed a dust 
sublimation temperature of 1700 K. All models with in- 
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Fig. 10. — Left: Reduced chi square map for the fit of our simple geometrical disk toy model (see text for description), as a function of 
two free parameters: inner disk radius and density power-law exponent. The flux ratio between the disk and the star has been optimized 
to minimize the \r values. Right: Optimum disk/star flux ratio associated to each disk model displayed on the x 2 map. 



put inner radii smaller than the sublimation radius have 
their inner radius forced to r- m = 0.1 AU and therefore 
give the same result. 

The small range of reduced \ 2 observed in the map of 
Fig. [TU] suggests that all models fit almost equally well 
our data set, so that the morphology of the circumstel- 
lar disk can not be meaningfully constrained. It must 
still be noted that models with very steep density pro- 
files (a > 6) starting at the sublimation radius have a 
slightly larger Xr- These particular models roughly cor- 
respond to narrow rings of dust located at the sublima- 
tion radius, which do not reproduce very well the data 
at the shortest baselines (see Fig. [!)]). For dust rings with 
larger inner radii, the oscillation frequency in the visi- 
bility curve is increased (the disk is fully resolved at a 
shorter baseline) and several local minima are then ob- 
tained in the x map, e.g. at r; n = 0.2 AU, 0.8 AU and 
1.5 AU. The solution at r; n = 0.2 AU would correspond 
to a sublimation temperature of 1300 K for black body 
dust grains. Conversely, models with flat density pro- 
files (a < 2) give equal fit qualities whatever the inner 
radius, because their spatially extended brightness dis- 
tribution produces almost no oscillation in the visibility 
curve (see green dash-dotted curve in Fig. [9]). It must 
also be noted that, whatever the chosen morphology of 
the dust disk, the best-fit disk/star flux ratio remains 
similar to the value derived with the uniform emission 
model. The largest discrepancy is found for the models 
with steep density power-laws starting at the sublima- 
tion radius, which give a best-fit contrast of 1.16% (i.e., 
about 2cr above our original estimation). This observa- 
tion reinforces our statement that the disk morphology 
has a weak influence on the best-fit disk/star contrast. 

6. DISCUSSION 

The discovery of hot dust grains in the inner plane- 
tary system of Fomalhaut cannot be regarded as a to- 
tally unexpected resul t. First, the data collected by 
iStapelfeldt et~aTl (pOOl ) with the MIPS and IRS instru- 
ments onboard Spitzer have already shown the presence 
of a warm excess emission at wavelengths between 17 
and 34 /im. MIPS imaging suggests that most of the 
24 yum excess is unresolved and originates from a 10" re- 
gion around Fomalhaut, i.e., well within the main dust 



ring located at 140 AU. Second, other A-type stars have 
already been shown to harbor similar near-infrared ex- 
cesses J^hjrtlravebejmin^ of hot 
dust (jAbsil et al.ll2006L 120081 : 1 Akeson et al.ll2009f ). In all 
previous cases, the absence of significant excess emission 
at wavelengths around 10 fim has suggested that the dust 
surface density distribution is very compact and does 
not follow a classical r -15 power-law. Modeling these 
interferometric observations has still shown that realis- 
tic debris disk models can be found to fit the observed 
near-infrared excesses of a few percent around A-type 
stars without producing too large amounts of excess at 
longer wavelengths, that would have been noticed in mid- 
infrared spectro-photometric measurements. 

It is not our purpose here to thoroughly study the prop- 
erties of the Fomalhaut inner disk, as this will be done in 
a companion paper, where we include additional interfer- 
ometric data collected at 10 (im with the Keck Interfer- 
ometer Nuller. It is nevertheless useful to check whether 
standard disk models could both reproduce the measured 
if-band excess and be consistent with photometric data 
at longer wavelengths. For that, we have used our debris 
disk toy model based on pure black body grains, assum- 
ing an r~ 15 surface density power law from the dust 
sublimation radius up to 140 AU. In order to reproduce 
the measured if-band excess, this disk must have a bolo- 
metric luminosity ratio idisk/^* ~ 5 x 10 -4 with respect 
to the stellar photosphere. The excesses produced at 17.5 
and 24 fim would then be of 0.38 and 0.22 Jy respectively, 
which would correspond to the measured IRS "on-star" 
excess of 0.4 ± 0.2 Jy at 17.5 /xm, while being fully within 
the me asured MIPS "on-star" excess of 0.6 ± 0.2 Jy at 
24 (Mm (jStapelfeldt et al.ll2004h . The predicted flux in 
the "extended" disk located outside 10", amounting to 
0.06 Jy, would also be within the MIPS PSF-subtracted 
estimation of 0.13 ± 0.02 Jy for the outer disk regions. 
Our model is also consistent with published photomet- 
ric measurements in the mid-infrared (including an IRAS 
12 /im observation), which are generally of poor accuracy. 
We will however show in a companion paper that this 
simple model cannot reproduce interferometric data ob- 
tained in the mid-infrared with the Keck interferometer, 
as well as unpublished IRS data in the 10-18 /im region. 
This suggests that the disk cannot be continuous under 
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the r" 1 ' 5 power law assumption, and therefore cannot 
be directly connected with the outer belt using classical 
models (e.g., based on Poynting- Robertson drag). 

7. CONCLUSION 

In this paper, we have re-analyzed published and un- 
published archival data from the VLTI/ VINCI instru- 
ment in search for near-infrared circumstellar emission 
around the young A4V star Fomalhaut. The large num- 
ber of visibility measurements at short baselines (146 in- 
dividual measurements with projected baselines smaller 
than 60 m), and the good accuracy of individual data 
points (median relative error bar of 2%), have allowed 
us to directly detect the presence of resolved circumstel- 
lar emission at a level of 0.88% ± 0.12% with respect 
to the stellar photosphere in the K band. If this ex- 
cess emission was produced by a dust disk similar to the 
solar zodiacal disk, its surface density would be about 
5000 times as large as in the solar system. Our analysis 
of the whole VINCI data set also provides an updated 
mean limb-darkened diameter 6*ld = 2.223 ± 0.022 mas 
for the stellar photosphere, which is predicted to have 
an apparent oblateness of 1.021. This measurement does 
not significant l y diffe r from the previous estimation by 
iDi Folco et alj (|2004f ). which did not take into account 
the presence of circumstellar emission. 

We attempted to constrain the morphology of the cir- 
cumstellar emission source with our interfcromctric data 
set, but failed to discriminate between a point-like source 
and an extended source, mostly due to the inappropriate 



observing strategy of this heterogeneous data set. Addi- 
tional constraints obtained with complementary observ- 
ing techniques have thus been used to determine the na- 
ture of the excess emission, showing that the presence 
of a point-like source within the interferometric field-of- 
view (either a bound companion or a background object) 
is very unlikely to reproduce the observed if-band ex- 
cess. Non-thermal extended sources are also shown to be 
mostly inconsistent with other observations. We there- 
fore consider that a hot dust population located mostly 
within the first 6 AU around Fomalhaut is the most likely 
explanation for the observed A'-band excess. 

In a companion paper, we will introduce further con- 
straints on the inner dust disk thanks to high dynamic 
range mid-infrared observations with the Keck Interfer- 
ometer Nuller, and complement these new constraints 
with various spectro-photometric measurements in order 
to investigate the nature and distribution of the dust 
grains. 
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APPENDIX 

OBLATE LIMB-DARKENED PHOTOSPHERIC MODEL 

To properly fit our interferometric data, which are taken at various baseline azimuths, we must take into account the 
influence of stell ar oblateness due to the rapid rotation of Fomalhaut (v sin i = 93 km s _1 , iRover et al.ll2007T ) . Following 
the discussion of I Absil et alj (|2008f ). we compute the apparent oblateness p of the Fomalhaut photosphere, defined as 
the ratio of the major and minor apparent radii (p = R a /Rb), using the following equation: 

where R m = (RgR b) 1 ^ 2 is the geometric mean between the major and minor apparent radii. Taking for R rn the stellar 
radius derived by IDi Folco et all (l2004h based on long-baseline interferometric measurements (R m = 1.840 i?©), and 
with M* = 2.0M ~(|Di Folco et al.l l2004). the apparent oblateness of the photosphere is p — 1.021. The orientation 
of the photosphere on the plane of the sky has recently b een constrained by spectrally dispersed interferometric 
measurements using VLTI/AMBER (Lc Bou auin et al.ll2009T ). showing a perfect m atch between the position angle of 
its major axis (PA star = 156°±3°) and that of the debris ring (PA disk = 156?0±0?3. lKalas et al.ll2005h . We will further 
assume that the photos pheric limb da rkening can be approximated by a linear law and use the i^T-band limb-darkening 
coefficient tabulated bv IClaretl (|2000h for an A4V star (uk = 0.19). With these assumptions, we are left with only 
one parameter to fit the interferometric data set: the mean limb-darkened angular diameter #ld- 

IMPORTANCE OF SHORT BASELINES FOR DEBRIS DISK DETECTION 

Another representation of the fit result described in Section 13.11 and Fig. [6] is proposed in Fig. [7J where we have 
separately displayed the data collected with the two shortest baseline configurations (E0-G0 and D0-H0). The squared 
visibility associated with the best-fit stellar photosphere model (fitted on the whole data set) is depicted by a blue 
curve, whose thickness represents the uncertainty on the squared visibility resulting from both the statistical and 
systematic error bars on its diameter. The presence of the circumstellar disk is revealed by the fact that all data 
points are located significantly below the expected stellar visibility. Even though the circumstellar disk is most evident 
at the shortest baseline (E0-G0), where the stellar photosphere is completely unresolved, the intermediate D0-H0 
baseline alone would be sufficient to conclude on its presence, although with a larger error bar on the final disk/star 
contrast. This figure clearly illustrates the importance of short baselines in the detection of circumstellar disks around 
main sequence stars, as the visibilities at longer baselines are much more affected by the uncertainties on the stellar 
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photo-spheric model. In particular, we stress that, thanks to the very short baselines used here, modifications in 
the fixed stellar limb-darkening and oblateness parameters would not significantly change the final estimation of the 
disk/star contrast. 
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